Luttinger liquid of polarons in one-dimensional boson-fermion mixtures 
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We use bosonization approach to investigate quantum phases in mixtures of bosonic and fermionic 
atoms confined in one dimensional optical lattices. The phase diagrams can be well understood in 
terms of polarons, which correspond to atoms that are "dressed" by screening clouds of the other 
atom species. For a mixture of single species of fermionic and bosonic atoms we find a charge 
density wave phase, a phase with fermion pairing, and a regime of phase separation. For a mixture 
of two species of fermionic atoms and one species of bosonic atoms we obtain spin and charge 
density wave phases, a Wigner crystal phase, singlet and triplet paired states of fermions, and a 
phase separation regime. Equivalence between the Luttinger liquid description of polarons and the 
canonical polaron transformation is established and the techniques to detect the resulting quantum 
phases are discussed. 



Mixtures of ultra-cold bosonic and fermionic atoms, 
that have recently become accessible experimentally, rep- 
resent a promising new system for studying strongly cor- 
related many-body physics [lj. Bosonic atoms mediate 
interactions between fermions and allow efficient cool- 
ing of the system Q. Several novel phenomena have 
been predicted theoretically for Boson-Fermion mixtures 
(BFM) including pairing of fermions formation of 
composite particles 0, spontaneous breaking of trans- 
lational symmetry in optical lattices |5| and appearance 
of charge density wave (CDW) 0. Most of these theo- 
retical studies relied on integrating out bosonic degrees 
of freedom to obtain an effective interaction between 
fermions, and then using a mean-field approach to in- 
vestigate many-body states 0- This approach, however, 
becomes unreliable in the regime of strong interactions. 
In particular, it fails in low-dimensional systems due to 
enhanced fluctuations and non-perturbative effects of in- 
teractions. 

In this paper we use bosonization method 0, |j| to 
investigate one dimensional (ID) BFM. The resulting 
quantum phases can be understood by introducing po- 
larons, i.e. atoms of one species surrounded by screen- 
ing clouds of the other species. Such dressed quasi- 
particles exhibit effective interactions and modified effec- 
tive masses. In our analysis the polarons emerge as quasi- 
particles with the slowest decaying correlation functions 
while quantum phases of the system arise from a com- 
petition of various ordering instabilities of such polarons. 
The phase diagrams we obtain (Figs. Ill-il) show a remark- 
able similarity to the Luttinger liquid phase diagrams of 
ID interacting electron systems |llj, suggesting that ID 
BFM may be understood as Luttinger liquids of polarons. 

In Fig. H we show a phase diagram for a mixture 
of bosons and spinless fermions as a function of exper- 
imentally controlled parameters: the scattering length 
between bosons and fermions (<x&/) and the strength of 
the longitudinal optical lattice for bosonic atoms 
|l2j . For relatively weak boson-fermion interactions and 
slow bosons (i.e. strong optical lattice for bosonic atoms) 
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FIG. 1: Phase diagram for a mixture of bosonic and spin- 
less fermionic atoms in a ID optical lattice. Shading in the 
/-PP phase describes the strength of the bosonic screening 
cloud (2A, see Eq.@) around a pair of fermions. Xl and Er 
are respectively the lattice period and recoil energy. Other 
parameters used for this figure are (see text for notations): 
v b =4,v f = 0.5, V b ,± = V f ,± = 20E R , V u = 2E R , boson- 
boson scattering length a bb = O.OIAl. 



the system is in the CDW phase, in which the densities 
of fermions and bosons have a 2fc/-modulation. In the 
case of very strong boson-fermion interactions the sys- 
tem is unstable to phase separation (PS) [l3l fl4j|. The 
two regimes are separated by a p-wave pairing phase of 
fermionic polarons (/-PP). While we carry out the de- 
tailed analysis for atoms in optical lattices, qualitatively 
it also applies to the continuous case 0, LJ liil • 

Before proceeding we note that bosonization approach 
has been applied to BFM in Ref. 0] . However this work 
did not consider correlations of polaronic degrees of free- 
dom, and as a result did not predict most of the quantum 
phases. The present system also resembles ID electron- 
phonon systems discussed previously in Ref. |is| . The key 
difference between the two is that the sound velocity in 
a solid state is typically smaller than the Fermi velocity, 
whereas for the dilute fermionic gas coupled to superfluid 
bosonic ensemble the opposite is true. In particular, this 
rules out an adiabatic elimination of the bosonic field. In 
contrast, we demonstrate in the present system polaronic 
quasi-particles dominate the low-temperature behaviour 
of the system. We also note that the 1-D p-wave super- 
fluid we obtain here may be of relevance to recent studies 
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in quantum information in Ref. Q . 

We first consider a mixture of spinless fermionic (/) 
and bosonic (6) atoms. For sufficiently strong optical 
potential the microscopic Hamiltonian is given by a single 
band Hubbard model 



(ij) * 
+ ~Y ^2 nb A n b,i - 1) + Ubf ^2 n b,i n f,ii (1) 



Here nb t i = b\b i and rifj = fjf i are respectively the 
boson and fermion density operators, and At&// are their 
chemical potentials. The tunneling amplitudes and 
the particle interactions Ub and Ubf depend on the laser 
beam intensities and the s-wave scattering lengths, and 
can be calculated explicitly (see e.g. Ref.[15j). In this 
paper we assume that the filling fraction of fermions 
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is not commensurate with the lattice or 



with the filling fraction of bosons iv The Fermi mo- 
mentum and velocity are given by kf = nisf and Vf = 
2tf sin(fc/), respectively. 

The Haldane's bosonization representations for 
fermion and boson operators are constructed as 0, Q 
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and 



m = 

b(x) = \v h + U b } 1/2 Y l ^- oc e 2mMb e l ' p >', where x is a 
continuous coordinate that replaces the site index i in 
Eq. (JJJ. The operators Hf/b(x) and $f/b(x) describe 
the bosonized density and phase operators for each 
of the atomic species, and satisfy the commutation 
relations \Ilb/f( x )> &b/f( x ')] — iS(x — x'). Short-ranged 
fluctuations are taken into account via the harmonics 
of Qf/b(x) = nvf/bx + ir J x dyILf/ b (y). The low energy 
effective Hamiltonian can be written to be the sum of 
the following terms: 
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(5) 



where we neglected the Umklapp scattering and back- 
ward scattering of bosons by assuming that bosons are 
deep inside the superfluid regime ^t|- When the boson- 
boson interaction is weak (i.e. 7 = m* h Ublv s < 10), 
the phonon velocity Vb and the scaling exponent of 
bosons Kb in Eq.Q can be well approximated by 
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where is the effective boson mass in the pres- 
ence of the lattice potential and v a ~ Vb is 



the superfluid fraction. The term in Eq. Q 

describes forward scattering between bosons and 
fermions, and the term H 2 ^ corresponds to back 
scattering of fermions on bosons (We defined fif(x) — 



c ft 

Jr. 



I R-p+kfh-p^ e 1X fL,p~kfR,p 

where ft R L j are the creation operators for the right 
and left moving fermions near the Fermi energy). 
It is useful, however, to take into account the de- 
viations of the phonon dispersion from the linear 
spectrum by using the Bogoliubov approximation 
w fc = y/(e b ,k - £b,o)(£b,k - £fc,o + 2UbV h ), where e b ,k is 
the band energy of a noninteracting boson. 

We integrate out the 2fc/-bosons and obtain an effec- 
tive fermion-fermion interaction. Since we are consider- 
ing a BFM with fast phonon velocity, i. e. ut > w/, we 
obtain within instantaneous approximation: 
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Here G 



92k f 



is the induced fermion-fermion interac- 



tion, and gk is the fermion- phonon (FP) coupling vertex, 
gk = Ubf \J Vb(£b,k — £&.o)/27rajfc. Note, that for small k 
we have a conventional FP coupling gk = glk] 1 ^ 2 with 
g = Ubf\fK~b~/'2ir. Therefore the effective Hamiltonian 
for a BFM is given by Eqs. ©-0 and @J with five pa- 
rameters: Vf, Vh, Kb, g and G. We can diagonalize this 
Hamiltonian [ijj and obtain 



dx 



j=a,A 



7Tn J (x) 2 + -(d x <i> J (x)y 

7T 



(7) 



where the eigenmode velocities, va and v a , are given by 
±(v 2 b +v})±y{v*-v}y + 16g*v b v f . (8) 
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Here Vf = (v 2 — AG) 1 / 2 and g = ge e with e e — ((vf — 
2G)/{v f + 2G)Y'\ 

When the FP coupling g becomes sufficiently strong, 
the eigenmode velocity va becomes soft, indicating an 
instability of the system to phase separation or collapse, 
depending on the sign of Ubf 0,0- To understand the 
nature of the many-body state of BFM outside of the 
instability region we analyze the long distance behavior 
of the correlation functions. For the bare bosonic and 
fermionic particles we find (b(x)b'(0)) ~ | | 2 an d 

(/(x)/t(0)) ~ cos(k f x)\x\-^ K ? +K y 1) (see Ref. ^0]). To 
describe particles dressed by the other species we intro- 
duce the composite operators 



fx{x) 



= e ~i\<S>b(x) 



f(x), b v (x)^e-^f^b(x), (9) 



with A and 77 being real numbers. The correlation func- 
tions of these operators are given by (f(x)f'(0)) ~ 
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cos(fc/a;)|xri( if ' 3+A2 ^ 1+ ^" 1 - 2A ^ 1 ) and (6(a:)&t(0)) ~ 
\ x \-U K 7 1 +'n 2K ^ 1 -^K-^ (geeRef.^3). We observe that 
the exponents of the correlation functions are maximized 
with A c = K e /K lt and r\ c = K 7 /K 7e , which we use to 
construct polaronic particles according to Eq. jjjj. In 
the limit of weak interactions we have A c — > Ubf /Ub and 
rj c — > 2Ub f I irvb- This result can be understood by a naive 
density counting argument that a fermionic polaron (/- 
polaron) locally suppresses a bosonic cloud by A c parti- 
cles, whereas a bosonic polaron (fe-polaron) depletes the 
fermionic system by rj c atoms. 

Ground states of one dimensional systems are often 
characterized by specifying the order parameters that 
have the slowest long distance decay of the correla- 
tion functions This is equivalent to finding the 
most divergent susceptibility in the low temperature 
limit (When the T — correlation function decays as 
(O(x)O(O)) ~ l/\x\ 2 - a , the finite T susceptibility di- 
verges as x{T) ~ l/T a ). For the 2k f CDW order param- 
eter, Ocdw — /l/r, we find acDW = 2 — 2Kp, and for 
the /-polaron pairing field, Of-pp = /lx/rx, we obtain 
a f -PP = 2-2 [X 2 K~ l + K~ l - 2XK~^\ . We did not in- 
clude polaron dressing in Ocdw, since this operator has 
no net fermionic charge and the exponent of Ocdw does 
not change if we replace / by /. For the /-PP operator, 
on the other hand, a/_pp is maximized when A = A c , 
so from now on ay_pp will always mean the /-PP expo- 
nent computed with the optimal A c . Scaling exponents 
shown in Fig. |2fa) demonstrate that divergencies of the 
CDW and /-PP susceptibilities are mutually exclusive 
and cover the entire phase diagram. In Fig. Efb) we 
show a global phase diagram of a BFM considering the 
FP coupling (g) and effective fermion-fermion interaction 
(G) as independent variables. This phase diagram is sim- 
ilar to what one finds for spinless electrons in Luttinger 
liquid theory The phase diagram in terms of ex- 
perimentally controlled parameters was shown in Fig. ^ 
We point out that if we were to define the pairing op- 
erator using bare fermions rather than /-polarons with 
optimal A c , we would have a region in the phase diagram 
in which none of the susceptibilities diverges (for the pa- 
rameters used in Fig. |2fa) this regime extends between 
0.55 < g/vf < 0.68). 

We next verify that the conventional construction 
of polaron operators based on the canonical polaron 
transformation (CPT) |2l| can give an equivalent po- 
laronic description as in Eq.[5]. The CPT operator is 
given by U\ = e ~ lX ^2k^o( Fkl3kP k +h - c -) j w here is the 
phonon annihilation operator, pk is the fermion den- 
sity operator, is some function of momentum k, 
and A specifies the strength of the phonon dressing. 
When applied to a fermion operator, the CPT gives 

U^f(x)U x = /(z)exp[-iA£ fc#0 (F fe /3 fe e- ife -*+h.c.) 
[2l| , which is the same as Eq. 10 , provided that one takes 




FIG. 2: Ground state of a BFM with spinless fermions. g is 
the longitudinal FP coupling and G is the effective fermion- 
fermion interaction after integrating out 2k f phonons. (a) 
Scaling exponents for Vb/vf = 3, Kb = 5 and G/t}f = 0.1. 
Different curves correspond to the 2k f CDW order parame- 
ter, /-polaron pairing field (/-PP), bare fermion pairing field 
(BFP), 6-polaron operator (b-P), and bare boson operator 
(BB). Note that operators constructed with polarons (bosons 
and Cooper pairs of fermions) always have larger exponents 
than their counterparts constructed with bare atoms, (b) 
Global phase diagram for Vb/vf — 5 and Kb = 10. 



Fk = |- v / ^ T sgn(fc) (note that in ID fermionic systems 
density operators correspond to Luttinger bosons). 

At finite temperature the correlation functions become 
(O(x)O(O)) - exp(-|x|/£)/M 2 ~ Q - The thermal correla- 
tion lengths for Ocdw and Of-pp are approximately 
given by £ ~ vj/kpT. For a finite system of length 
L = NXp (N being the number of lattice sites), the 
T = 0-properties of the system are visible for £ ~ L. 
This corresponds to a temperature regime of T ~ Tf/N, 
Tf being the Fermi temperature. 

Before concluding our discussion of BFM with spin- 
less fermions we discuss an approach for observing the 
phase transition between the CDW and the /-PP phases 
in laser stirring experiments |22j. Suppose that a laser 
beam is focused at the center of the cloud, such that it 
creates a weak local potential for fermionic atoms. In the 
/-PP phase the stirring potential can be moved through 
the system with no dissi pat ion, if its velocity is slower 
than some critical value 22]. At the /-PP/CDW phase 
boundary the critical velocity goes to zero, reflecting a 
transition to the insulating (CDW) state. This scenario 
follows from an RG analysis of a small local potential 
for fermions: Such a potential is irrelevant in the /-PP 
phase, but becomes relevant in the CDW phase, similar 
to the single impurity problem of a ID electron system 
|23j |. It should be emphasized that only when polarons 
are used to construct Cooper pairs, the paired phase co- 
incides with the domain of irrelevance of a small pinning 
potential. 

We now extend our analysis to BFM with fermions 
with two internal hypcrfine states, which we assume 
to be SU (2) symmetric. Spin symmetry of the sys- 
tem leads to separation of the bosonized Hamiltonian 
into spin and charge sectors. The charge part of 
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FIG. 3: Phase diagrams for a mixture of bosonic and S — 1/2 
fermionic atoms. Shading in the paired regions describes the 
strength of bosonic screening clouds around fermionic po- 
larons. In both figures Vb/vj = 5 and Kb = 10. In (a) 
Ufi/vf = -0.8tt and in (b) U n /v f = 0.8tt. At the top of (b) 
and everywhere in (a) we have gix < (K a < 1), and gi± < 
(K a > 1) at the bottom of (b). For larger G and Ufi one can 
also have a Wigner crystal phase (not shown here) . Note that 
parentheses (. . .) indicate subdominant phases, while square 
brackets [. . .] indicate degenerate phases. 



the Hamiltonian is equivalent to a BFM with spinless 
atoms and can be diagonalized analogously. The spin 
part of the Hamiltonian has a form of the sine-Gordon 

model, H a = \v a J dx nU a (x) 2 + ± (d x $ a (x)) 

2gi 



(27^)2 ' / dxcos [y/8K a Q a (x)\ , where v a is the spin ve- 
locity, K a is the spin Luttinger exponent, and — 
U-\i — 4ttG is the effective backward scattering ampli- 
tude for fermions, that has contributions from the bare 
fermion-fermion interaction and from integrating out 2k f 
phonons. The nature of spin excitations in the ground 
state follows from the well known properties of the sine- 
Gordon model: For K a > 1 the system has gapless spin 
excitations (<?i_l is irrelevant), for K a < 1 the system has 
a spin gap (g\±_ becomes relevant). In order to describe 
possible ground states of the system we calculated the 
low temperature behavior of susceptibilities for the fol- 
lowing order parameters: 2k f spin density wave (xsdw), 
2kf charge density wave (xcdw), 4Jtf (Wigner crys- 
tal) charge density wave (xwc), singlet polaron pairing 
(xspp), and triplet polaron pairing (xtpp)- Depend- 
ing on the parameters we find the following regimes: 1) 
Equally divergent xcdw and xsdw- Degenerate CDW 
and SDW phases; 2) Equally divergent xspp and xtpp- 
Degenerate SPP and TPP phases; 3) Divergent XSPP- 
State with singlet pairing of polarons; 4) Both Xcdw 
and xspp diverge, but unequally. This probably corre- 
sponds to a phase that has both a CDW order and singlet 
polaron pairing; 5) Divergent xcdw- CDW phase; 6) Di- 
vergent xwc- Wigner crystal phase. The phase diagram 
in Fig. |3 again shows remarkable similarity to a phase 
diagram for interacting electrons [Tl| . 

In summary, we used the bosonization method to in- 
vestigate ID mixtures of bosonic and fermionic atoms 
involving spinless and S = 1/2 fermions. Interactions 
between atoms can lead to such interesting phenomena 



as spin and charge density waves, singlet and triplet pair- 
ing of atoms. This phase diagram can be understood in 
terms of polarons, in that it corresponds to a Luttinger 
liquid of polarons. We also discussed laser stirring exper- 
iments as a technique for probing quantum phase transi- 
tions between paired and insulating phases of polarons, 
and considered the finite temperature and size of a real 
system. 
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